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This paper presents a comparison of two finite-temperature BCS-Bose Einstein condensation (BEC) crossover 
theories above the transition temperature: Nozieres Schmitt-Rink (NSR) theory and finite T-extended BCS- 
Leggett theory. The comparison is cast in the form of numerical studies of the behavior of the fermionic spectral 
function both theoretically and as constrained by (primarily) radio frequency (RF) experiments. Both theories 
include pair fluctuations and exhibit pseudogap effects, although the nature of this pseudogap is very different. 
The pseudogap in finite T-extended BCS-Leggett theory is found to follow a BCS-like dispersion which, in turn, 
is associated with a broadened BCS-like self energy, rather more similar to what is observed in high temperature 
superconductors (albeit, for a d-wave case). The fermionic quasi-particle dispersion is different in NSR theory 
and the damping is considerably larger We argue that the two theories are appropriate in different temperature 
regimes with the BCS-Leggett approach more suitable nearer to condensation. There should, in effect, be little 
difference at higher T as the pseudogap becomes weaker and where the simplifying approximations used in the 
BCS-Leggett approach break down. On the basis of momentum-integrated radio frequency studies of unpolar- 
ized gases, it would be difficult to distinguish which theory is the better A full comparison for polarized gases 
is not possible since there is claimed to be inconsistencies in the NSR approach (not found in the BCS-Leggett 
scheme). Future experiments along the lines of momentum resolved experiments look to be very promising in 
distinguishing the two theories. 
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I. INTRODUCTION 

The behavior of ultracold superfluid Fermi gases contin- 
ues to attract attention from the experimental and theoretical 
communities. Through a Feshbach resonance iH 01 which 
tunes the strength of the attractive interaction, these trapped 
gases can exhibit a crossover from BCS to Bose Einstein con- 
densation (BEC). Unlike the Bose superfluids, at the time of 
their discovery, the Fermi superfluids are not associated with 
any ready-made theory (such as Gross Pitaevskii, or Bogoli- 
ubov theory for bosons). This provides an opportunity for 
theorists to work hand in hand with experimentalists to arrive 
at the same level of understanding of the fermionic as was 
reached for the bosonic superfluids. While the Fermi systems 
are harder to address than the Bose counterparts, the payoff 
for progress is great. Moreover, there is a general belief that 
these systems may lead to important insights into high tem- 
perature superconductor (HTSCs), in part because the HTSCs 
exhibit an anomalously short coherence length which suggests 
that they may be mid-way between BCS and BEC iHSbSI- 

There is, as yet, no clear consensus about the theory 
which underlies BCS-BEC crossover although there are three 
rather well studied analytic many body theories which have 
emerged. The goal of the present paper is to present a com- 
parison and assessment of these approaches with a particular 
focus on the two of the three, which seem most reliable. In 
addition to assessing the theoretical approaches, we address 
physical consequences and how the theories may be differen- 
tiated through the behavior of the centrally important spectral 
function and related density of states. We do this in the con- 
text of radio-frequency-(RF) based probes. 

What is among the most interesting features of BCS-BEC 



crossover is the fact that the normal state (out of which super- 
fluidity forms) is different from the normal state (Fermi liquid) 
associated with strict BCS theory. The normal state of, for ex- 
ample, a unitary gas consists of pre-formed pairs which persist 
below Tc in the form of non-condensed pair excitations of the 
condensate. This excitation branch is in addition to the usual 
gapped fermionic excitations. The normal state is often said to 
exhibit a "pseudogap" which has features in common with the 
exotic normal state of the high temperature superconductors. 
This pseudogap [0, 01 reflects the formation of quasi-bound 
pairs which in turn require an energy input (called A) in order 
to break the pairs and create fermionic excitations. Physically, 
what differs from one crossover theory to another ||4j, j7ll8l|9|l is 
the nature of these non-condensed or pre-formed pairs which, 
respectively, appear below and above Tc- Unlike the pair fluc- 
tuations of traditional superconductors (which are associated 
with low dimensionality and impurity effects) these pairs are 
present because of a stronger-than-BCS attractive interaction. 
As a consequence, the pairing gap A persists to temperatures 
which can be several times Tc for the case of the unitary gases. 

In this paper we address the temperature dependence of the 
spectral function particularly in the normal state. The den- 
sity of states (DOS), which can be obtained from the spec- 
tral function, will also be presented. We compare with ex- 
periments in the context of RF spectra of both unpolarized 
and polarized Fermi gases. Quantum Monte Carlo simula- 
tions 1 10. 11] provide useful information such as the super- 
fluid transition temperature Tc, entropy, condensate fraction, 
etc. and recently reveal evidence of non-condensed pairs fl^ 
along with a pseudogap |13] in the normal phase. Our fo- 
cus is on two different finite-temperature BCS-BEC crossover 
theories and we present a detailed comparison of the results 
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obtained from the two theories as well as an assessment of 
other BCS-BEC crossover theories. 



1. Analysis of Different Crossover Theories 

A fair amount of controversy [0, 13 [lH [T^ has surfaced 
in the literature regarding the three alternative analytic pairing 
fluctuation schemes. In this paper we address some of these 
issues and clarify misleading claims. At this early stage of un- 
derstanding we do not believe it is suitable to invoke (possibly 
fortuitous) fits to particular experimental or Monte-Carlo de- 
rived numbers to establish which of these theories is "best". 
Rather in line with the goal of this paper, one has to look at 
the differences at a more general level. One has, furthermore, 
to subject these theories to careful consistency tests. 

Each of the three many body approaches is associated with 
a different ground state. Thus far, only one of these can be 
written down analytically. In this context we note that one 
can trace the historical origin of the BCS-BEC Uterature to 
the observation that the simplest BCS-like wavefunction 

= nk(wk + fk4^cLk|)|o}, (1) 

is much more general than originally presumed UtI [Tsl [l9ll . 
To implement this generalization, all that is required is that 
one solve for the two variational parameters Wk and Uk in con- 
cert with a self consistent condition on the fermionic chem- 
ical potential /i. As the attraction is increased, ^ becomes 
different from the Fermi energy Ep, and in particular, in the 
BEC regime, /i is negative. This ground state is often called 
the "BCS-Leggett" state and the two variational parameters 
Mk and Vk can be converted to two more physically accessi- 
ble parameters associated with the zero temperature gap (or 
equivalently order parameter) and ji. 

The three theories currently of interest can be related to a 
t-matrix scheme. Within a given t-matrix scheme one treats 
the fermionic self energy and the pair-propagator using a cou- 
pled equation scheme, but drops the contributions from higher 
order Green's functions. This t-matrix is called tpg{Q), where 
Q = {ifli, q) is a four- vector and fti denotes the boson Mat- 
subara frequency; it characterizes the non-condensed pairs 
which are described physically and formally in different ways 
in the different theories. Here the subscript pg is associated 
with the pseudogap (pg) whose presence is dependent on the 
non-condensed or pre-formed pairs. Quite generally we can 
write the t-matrix in a ladder-diagram series as 

U 

" i + UxiQY 

where xiQ) the pair susceptibility and U denotes the attrac- 
tive coupling constant. 

The Nozieres-Schmitt-Rink (NSR) theory 1 19] is associated 
with a pair susceptibility x{Q) which is a product of two bare 
Green's functions. The fluctuation exchange or FLEX ap- 
proach is associated with two dressed Green's functions and 
has been discussed by Haussmann ll20ll2l|] . Zwerger and their 
collaborators in the context of the cold gases, and even earlier 



in the context of the cuprates | [22l l23l l24ll . It is also called the 
Luttinger-Ward formalism lll6ll . or Galitskii-Feynman theory 
JH. FinaUy, it is well known llH IH that BCS theory 
(and now its BCS-BEC generalization) is associated with one 
bare and one dressed Green's function in the pair susceptibil- 
ity. 

These differences would seem to be rather innocuous and 
technical but they have led to significant qualitative differ- 
ences and concurrently strong claims by various proponents. 
We stress that while there are several variants, as we dis- 
cuss below, the version of the NSR scheme which seems to 
us most free of concerns is that discussed in References |j7|] 
which introduced a more physical treatment of the number 
equation. This revision of strict NSR theory was, in part, an 
answer to J. Serene 129(1 who raised a question about a cen- 
tral approximation in the theory in which the number equation 
(n = 2 G(K), where G{K) is the single particle Green's 
function) is approximated by n = — ^g^'' , where the thermo- 
dynamical potential Vtth is approximated by a ladder-diagram 
series. It was shown that this amounts to taking the leading 
order in a Dyson series for G{K). 

The present paper concentrates on the normal state behav- 
ior, although all three classes of theories have been extended 
below Tc- What is essential about these extensions is that the 
non-condensed pair excitations associated with tpg are gap- 
less, as in boson theories. Indeed, it is in these T <Tc exten- 
sions that a number of concerns have been raised. In partic- 
ular, in the leading order extended NSR theory (or so called 
"Bogoliubov level" approach), ll9l [30l fsill . the gap equation 
(which is assumed to take the usual BCS form, rather than 
derived, for example, variationally) does not contain explicit 
pairing fluctuation contributions; these enter indirectly only 
via the fermion chemical potential /i. At this level, the number 
equation is the only way in which explicit pairing fluctuations 
are incorporated. At the so called "Popov level" calculation, 
the gap equation is presumed to contain pair fluctuations [3j] 
but there is some complexity in ensuring the concomitant gap- 
lessness of the pair excitations. Similar issues arise with the 
FLEX or Luttinger-Ward approach in which (| 16] and refer- 
ences therein) gapless sound modes must be imposed some- 
what artificially. 

While the order of the transition at Tc is second order in 
the BCS-Leggett scheme it is first order ]11,|1] in NSR based 
approaches (as well as for the fully renormalized pair suscep- 
tibility scheme. This leads to unwanted features in the density 
profiles I's^ and T dependent superfluid density ]33], Ps{T)- 
Despite these unphysical aspects, the NSR-based scheme cap- 
tures the physics of Bogoliubov theory of weakly interact- 
ing bosons and should, in principle, be the quantita- 
tively better low T state, particularly in the BEC limit. Nev- 
ertheless some issues have been identified ll34ll which sug- 
gest the breakdown of true quasi-particles associated with 
Bogoliubov-like theories for paired fermions. This, in turn, 
derives from the self consistent treatment of coupling between 
the non-condensed pairs and the sound modes. Further anal- 
ysis will be required to establish if this is compatible with 
experimental or theoretical constraints. 

A very early concern about the so-called "GG" or FLEX 
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approach was raised in a paper by Kadanoff and Martin |126 
in 1961: "The similarity [to a Bethe-Salpeter equation] has 
led several people to surmise that the symmetrical equation 
[involving fully dressed G's everywhere] solved in the same 
approximation would be more accurate. This surmise is not 
correct. The Green's functions resulting from that equation 
can be rejected in favor of those used by BCS by means of a 
variational principle." Importantly this approach does not have 
a true pseudogap. Despite claims by the Zwerger group ITgIi 
that theirs is a more fully "consistent" theory, and in this con- 
text appealing to Ref. Il35ll . the authors of Ref. |35] instead 
say: "We thus conclude that ... approaches such as FLEX 
are unreliable in the absence of a Migdal theorem and that 
there is indeed a pseudogap." Similar observations have ap- 
peared elsewhere in the literature As noted 
in Ref. [36] " vertex corrections to the self energy, which are 
discarded in the previous studies [of FLEX] are crucially im- 
portant for the pseudogap". Additional concerns have been 
noted recently [25] that in the FLEX (or GG t-matrix) theory 
the propagator G does not display quasiparticle poles asso- 
ciated with the gap. "This is because the Dyson equation, 
G{k) = l/(z — k^/2m — excludes identical poles of 
G and E while the linear relation demands them". 

In recent work below |^ [sTl Ultj a non-variational gap 
equation was used to derive an additional term in the number 
equation related to dVtth / d/S.sc 0. Here A^c is the order 
parameter and, here, again, Q,th is the thermodynamical po- 
tential. This extra term means there is no variational free en- 
ergy functional, such as required by Landau-Ginsburg theory. 
Of concern are arguments that by including dVlth / dAgc ^ 0, 
it is possible to capture the results of Petrov et al [38] for 
the inter-boson scattering length. We see no physical con- 
nection between the exact four-fermion calculations and the 
non-variational component of the many body gap equation. It 
should, moreover, be stressed that all other t-matrix schemes 
have reported an effective pair-pair scattering length given by 
ub — 2a which is larger than the value ub = 0.6a obtained 
from a four-body problem Ii38i] . Here a is the s-wave scatter- 
ing length of fermions. Indeed, our past work ll39ll and that 
of Reference |32] have shown that one needs to go beyond 
the simple t-matrix theory to accommodate these four-fermion 
processes. 

Additional concerns arise from the fact that an NSR-based 
scheme has difficulty I40„ i41il accommodating polarization ef- 
fects in the unitary regime. As stated by the co-workers in 
Reference (^Tj: "Unfortunately, in a region around the uni- 
tary limit we find that the NSR approach generally leads to a 
negative population imbalance at a positive chemical potential 
difference implying an unphysical compressibility.". 

The central weakness of the BCS-Leggett approach (and 
its finite-T extension) appears to be the fact it focuses prin- 
cipally on the pairing channel and is not readily able to 
incorporate Hartree effects. The evident simplicity of this 
ground state has raised concern as well. Clearly, this is 
by no means the only ground state to consider but, among 
all alternatives, it has been the most widely applied by the 
cold gas community includin g th e following notable papers 
llil H El m Si m, Mm MMM- The central 



strengths of the finite-T extended BCS-Leggett approach in 
comparison with others are that (i) there are no spurious first 
order transitions and (ii) the entire range of temperatures is 
accessible, (iii) Moreover, polarization effects may be read- 
ily included 1.53. ,54.1 , (iv) as may inhomogeneities which are 
generally treated using Bogoliubov deGennes theory lIsoifsTll . 
based on this ground state. 

The above analysis leaves us with two theoretical schemes 
which we wish to further explore: the NSR approach (which 
in the normal phase follows directly from the original paper 
lfl9ll ) and the BCS-Leggett-based scheme, as extended away 
from zero temperature, and in particular above Tc- As t-matrix 
approaches to many body theory, these are similar in spirit, 
but different in implementation. It is clearest below Tc, that 
the two theories focus on different physics. NSR approaches 
view the dominant processes as the coupling of the order pa- 
rameter collective modes to the non-condensed pairs and the 
BCS-Leggett scheme focuses on the steady state equilibrium 
between the gapped fermions and the non-condensed pairs. 
Thus NSR focuses more fully on the bosonic degrees of free- 
dom and BCS-Leggett focuses on the fermionic degrees of 
freedom. Above Tc, because the NSR scheme involves only 
bare Green's functions, it is simpler. Thus, it has been studied 
at a numerical level in a more systematic fashion. In the litera- 
ture, the BCS-Leggett a ppro ach at T ^ 0, has been addressed 
numerically ll25ll55l l56', 57], assessed more theoretically |22|, 
as well as applied to different physical contexts [58, .59„ ,60ll . 
In this paper we apply an approximation based on prior nu- 
merical work ll55il56ll to simplify the calculations. 



2. The Fermionic Spectral Function 

A central way of characterizing these different BCS-BEC 
crossover theories is through the behavior of the fermionic 
spectral function, A(k, uj). For the most part, here, we restrict 
our consideration to the normal state where A(k, oj) should 
indicate the presence or not of a pseudogap. A momentum 
integrated form of the spectral function is reflected in radio 
frequency studies- both tomographic |6lll or effectively ho- 
mogeneous and trap averaged 11621 16311 . One of the princi- 
pal observations of this paper is that these momentum inte- 
grated probes are not, in general, sufficiently sensitive to pick 
up more than gross differences between the three crossover 
theories. 

However, there are now momentum resolved RF studies 
ll64|l which probe the spectral function more directly, in a 
fashion similar to angle-resolved photoemission spectroscopy 
(ARPES) probes of condensed matter A central aim of this 
paper is to show how these studies in future will be able 
to differentiate more clearly between the different crossover 
schools. Here we confine our attention to homogeneous 
systems, although experiments are necessarily done for the 
trapped case. In addition to RF spectroscopy, it was proposed 
[65] that the spectral function can also be measured in Raman 
spectroscopy. 

We note that for the HTSCs, ARPES studies have been cen- 
trally important in revealing information about the supercon- 
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ducting as well as the pseudogap phases ||66|] . Indeed, the 
close relation between ARPES and radio frequency probes 
has been discussed in our recent review ll67ll . It was shown 
in Ref. (6^ that the spectral function of HTSCs in the pseu- 
dogap phase appears to exhibit dispersion features similar to 
those in the superconducting phase. This spectral function is 
modeled ifssl O^L l69ll by a broadened BCS form with a self 
energy 



Ap,(k)2 



(2) 



Here Apg(k) is the (s or d-wave) excitation gap of the nor- 
mal phase and 7 is a phenomenological damping. Frequently, 
one adds an additional, structureless imaginary damping term 
iEo, as well. High temperature superconductor experiments 
at temperatures as high as T « l.STc have reported that in the 
regions of the Brillouin zone (where the pseudogap is well es- 
tablished), the dispersion of the fermionic excitations behaves 
like 



i^k«±J(ek-Ai)' + Ap<,(k)2 



(3) 



Importantly Eq. ^ has also been used in the cold gas stud- 
ies 1 64] in the region near and above Tc and implemented phe- 
nomenologically below Tc ll65ll . This both demonstrates the 
presence of pairing and ultimately provides information about 
the size of the pairing gap. It has been shown that Eq. (|2]i 
is reasonably robust in the extended BCS-Leggett state above 
Tc, at least up to temperatures |'55','56'l of the order of w l.STc- 
By contrast this approximate self energy is not generally suit- 
able to NSR theory |55, 56], although for T/Tc = 1.001 a fit 
to Eq. (|3]l has been obtained. In a similar context we note that 
in the FLEX approach, the spectral function and associated 
self energy is not of the broadened BCS form. Mathemati- 
cally, this BCS-like structure in the self energy and fermionic 
dispersion (which is numerically obtained) comes from the 
facts that the effective pair chemical potential ^pair — > at 
and below Tc, and that by having one bare and one dressed 
Green's function in xiQ) there is a gap in the pair excitation 
spectrum so that the pairs are long lived; in this way 7 (which 
scales with the inverse pair lifetime) is small. Physically, we 
can say that this behavior reflects the stability of low momen- 
tum pairs near Tc and below. 

These differences between the three different crossover the- 
ories become less apparent for the momentum integrated RF 
signals. In the BCS-Leggett approach at low temperatures 
the dominant structure comes from pair breaking of the con- 
densate (which would be associated with the negative root 
in Eq. (O). Despite the fact that their fermionic dispersions 
are different, both other theories yield a very similar "pos- 
itive detuning branch" in the RF spectrum llm FtoIi . How- 
ever, at higher temperatures, both for polarized and unpolar- 
ized gases, there is theoretical evidence of the "negative de- 
tuning branch" arising from the positive root in Eq. (O in the 
BCS-Leggett based approach ll67l FtTIi . This is absent in the 
two other schemes, at least within the normal state. It also ap- 
pears to be difficult to see experimentally in the unpolarized 



case, although it is clearly evident once even a small polariza- 
tion lis m is present 172]. 

This paper is organized as follows. Sections |ll] and |III] 
briefly review NSR theory and BCS-Leggett theory as ex- 
tended to non-zero T. Section |IV] addresses a comparison of 
the spectral function at unitarity and on the BEC side obtained 
from the two theories. In subsections |IV B| we plot a compar- 
ison of the related density of states at unitarity and in IIV CI 
we address a comparison of RF spectra in the two theories for 
an unpolarized Fermi gas which also addresses experimen- 
tal data. Also included is a prediction of RF spectra on the 
BEC side of resonance. The remaining sections (Section [Vl 
and Section rvib do not focus on comparisons because the is- 
sues discussed pertain to questions which only BCS-Leggett 
theory has been able to address. Here we propose a subtle sig- 
nature of the superfluid transition in Section [V] which could 
be addressed in future and in Section[yT]we address the theo- 
retical RF spectrum of polarized Fermi gases at unitarity and 
its comparison with experimental data. Section lVlIl concludes 
this paper. 

We remark that in this paper we study s-wave pairing in 
three spatial dimensions which is more relevant to ultra-cold 
Fermi gases while HTSCs should be modeled as d-wave pair- 
ing in quasi-two dimensions. However, as one will see, there 
are many interesting common features in these two systems. 



n. NSR THEORY ABOVE Tc 

The normal state treatment of NSR theory which we ap- 
ply here follows directly from the original paper in Ref. 111911 . 
Here the different variants of NSR theory as introduced by 
different groups ^ [isl [13, [32|l are not as important. Al- 
though there is still the concern 129] that the number equation 
is only approximate, the numerics are simpler if we follow 
the original approach; comparisons with more recent work 
in Ref. 173] (based on the fully consistent number equation 
(n = 2 J2k C!{K)) seem to validate this simplification. This 
same more consistent number equation is used throughout the 
work by the Camerino group |0] • 

NSR theory builds on the fact that the fermion-fermion at- 
traction introduces a correction to the thermodynamic poten- 
tial: A^th = flth-^f = I]Qln[?7"^+xo(Q)], where ri/ is 
the thermodynamic potential of a non-interacting Fermi gas, 
XoiQ) is the NSR pair susceptibility, and J^q — Tj^i Sq- 
In the normal phase. 



Xo 



(0) = J2Go{K)GoiQ - K) 



K 



E 



/(ek+q/2) + /(a-q/2) - 1 

i^l - (Ck+q/2 + Ck-q/2) 



(4) 



Here K = {ioj^, k) where lOi, is the fermion Matsubara fre- 
quency, Y.K = ^E,.Ek' Go(i^) = l/{iuju - Ck) is the 
non-interacting fermion Green's function, = k"^ /"^m — ji 
where m and fi denote the mass and chemical potential of the 
fermion, and f{x) is the Fermi distribution function. We set 
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tpg 

AAAAAV = 



\ (T^Tc) 



Figure 1 : Fermion self-energy and ladder diagrams in (a) NSR theory 
and (b) GGo t-matrix theory. The thin solid line, thick solid line, 
thick dashed line, dotted line, wavy line, vertical thin dashed line 
denote Go, G, to, tsc = -AL(5(Q)/r, tpg, U. 



Here toR{^, q) = to{i^i ^ + iO"'"; q) and V denotes the 
Cauchy principal integral. In the presence of bound states, 
tQii{n, q) has poles which result in bound state contributions 
to the fermion self-energy 



q 



dRelt- 
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9Re[t-^] 









1 



Qbs — UJ — Cq-k_ 
6{Qbs - UJ - ^q-k) 



(7) 



Here ilbs = fibs (q) denotes the location of the pole in toi?- 



ni. BCS-LEGGETT THEORY: BROKEN SYMMETRY 
PHASE 



h = 1 and fcs = 1. NSR theory is constrained by the con- 
dition U^^ + xo(0) > since U'^ + xo(0) = signals 
an instability of the normal phase and the system becomes a 
superfluid as temperature decreases. 

The fermion chemical potential is determined via the NSR 
number equation 



(5) 



where n is the total fermion density. As noted in fW], when 
/i < 0, +xo{i^i ri + iO+,q) = Ohas solutions which 
correspond to bound states. Those bound states have contri- 
butions proportional to b{flq'') to the total density, where b{x) 
is the Bose distribution function and ilq"-* is the energy disper- 
sion of the bound states. 

The NSR Green's function is G(K) = [Go{K)-T.o{K)]-'^ 
and its retarded form is Gii{u!,'k.) = G{iujn ^ lo + 
«0"'",k). Following a t-matrix formalism, one can also con- 
sider the coiTections to the fermion self-energy So(i^) = 
X]q ^o(Q)Go(Q — K). Figure [T] illustrates the structure of 
fermion self-energy in NSR theory and in the finite temper- 
ature theory associated with the BCS-Leggett ground state, 
which will be summarized in the next section. Here the t- 
matrix is given by to(Q) = ^IW^^ + Xo(Q)]- The retarded 
form of the fermion self-energy has the structure Eo(jw,/ — > 
w + iO+, k) = E;,(cj, k) + iT.'^{uj, k), where Sf, and Sf,' cor- 
respond to the real and imaginary part of the self-energy. We 
separate the contribution of the bound states from the rest 
(called the continuum contribution). The continuum contri- 
bution is 

^',A^,Lo) = ^{/(Cq)Re[ioii(q + k,c. + Cq)] + 



V 



dn b{n) 



-Im[tofl(q + k,r!)]|, 



We first review BCS-Leggett theory as it has been applied in 
the broken symmetry phase. The first three equations below 
represent a t-matrix approach to the derivation of the stan- 
dard BCS gap equation. In this way we set up a machinery 
which is readily generalized to include BCS -B EC crossover 
theory. BCS theory can be viewed as incorporating virtual 
non-condensed pairs. Because they are in equilibrium with 
the condensate, the non-condensed pairs must have a vanish- 
ing "pair chemical potential", fXpair = 0. Stated alternatively 
they must be gapless. The t-matrix can be derived from the 
ladder diagrams in the particle-particle channel (see Fig.[T]i: 



tpgiQ) 



u 



1 + UJ:kG{K)Go{-K + Q)' 



(8) 



with tpg{Q = 0) ^ oo, which is equivalent to Upair = 0, for 
T < Tc- Here G, and Gq represent dressed and bare Green's 
functions, respectively. To be consistent with the BCS ground 
state of Eq. ([T]), the self energy is 

= Y.^sc{Q)Go{-K + Q) 



Q 



Q 

-AlGoi-K). 



(9) 



Asc{T) is the order parameter while A(T) is the pairing 
gap. From this one can write down the full Green's func- 
tion, G{K) = [G^^{K) - i:sciK)]-\ Finally, Eq. dD with 
Mpair — gives the BCS gap equation below T^. 



1 



k 



1 " vjEr) 

2Eff 



(10) 



with = 



(6) 



+ Al^. We have, thus, used Eq. (O 
to derive the standard BCS gap equation within a t-matrix lan- 
guage and the result appears in Eq. ( fTOl ). Eq. ([8]) above can 
be viewed as representing an extended version of the Thou- 
less criterion of strict BCS which applies for all T < Tc- 
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This derivation leads us to reaffirm the well known result 
ll26l IztI I28I1 that BCS theory is associated with one bare and 
one dressed Green's function in the pair susceptibility. 

Next, to address BCS-BEC crossover, we feed back the con- 



tribution of the non-condensed pairs which are no longer vir- 
tual as they are in strict BCS theory, above. Eq. ([8]l is taken as 
a starting point. Equation (|9]l is revised to accommodate this 
feedback. Throughout, Q denote four-vectors. 



J 



T.{K) = tiQ)Goi-K + Q)= Y,[ts 
Q Q 



(Q) + tpg{Q)]Go{-K + Q) = Esc(K) + ^pg{K) 



(11) 



Numerically, Spg(if) 



A2 

pg 



iuoy + Ck - M + *7 



zEo; analytically, Yisc{K) 



A 



(12) 



7, So small: E(i^) « -(A^ + ^lg)Go{-K) ^ -A^Goi-K) A^^ ^ -^tp,(Q) 

Q 



(13) 



tpgiQ = 0) 



2/(£^k) 



I — 



Ek = v/(ek - + A2 



(14) 



Note that Eqs. (|9]l and (fTTT i introduce the self energy which is 
incorporated into the fully dressed Green's function G{K), 
appearing in tpg. Also note the number equation n = 
2 J2k G{-^) is solved consistently: 



K 



^k -t^k 



(15) 



where ^k = Ck ^ 

This leads to a closed set of equations for the pairing gap 
A(r), and the pseudogap Apg(r) (which can be derived 
from Eq. (fT3]l). The BCS-Leggett approach with the disper- 
sion shown in Eq. ( fT4b thus provides a microscopic deriva- 
tion for the pseudogap model implemented in Ref. ll65ll . To 
evaluate A^g (T) numerically, we assume the main contribu- 
tion to tpg{Q) is from non-condensed pairs with small Q, 
which is reasonable if temperature is not too high ifssl ISol . 
The inverse of tpg after analytical continuation is approxi- 
mated as tpg{uj, q) « [Z{il — ilq + Upair) + "i^q]^^^ where 
Z = (9x/5r!)|j2=o,g=o, = qV{2Mb) with the effec- 
tive pair mass M^^^ = {l/3Z){d^x/dq'^)\ 0=0, g=o which 
takes account of the effect of pair-pair interactions. Near Tc, 
Tq — > faster than as q ^ 0. Following this approxima- 
tion, Apg{T) essentially vanishes in the ground state where 
A = A,c. 

The entire derivation contains one simplifying (but not fun- 
damental) approximation. Detailed numerical calculations 
llssi I56I1 show that T,pg can be written as in Eqs. (|2]l, which 
is the same as that in Eq. ( fT2] i. with the observation that as 
Tc is approached from above, 7 and Sq which appears in 
Eq. (fT2] i become small. To zeroth order, then, we drop 7 and 
So (as in Eq. (fTJIi), and thereby can more readily solve the 
gap equations. To first order we include this lifetime effect as 



in Eq. (fT2l i in addressing spectral functions and other correla- 
tions. 

The actual value of 7 makes very httle qualitative differ- 
ence and the previous numerical calculations ifssl l56ll do not 
include d-wave or trap effects so that we should view 7 as a 
phenomenological parameter. For the HTSCs, the expression 
for Tipg in Eq.© is standard in the field ll69[[74l p75ll. and we 
can use specific heat jumps or angle resolved photoemission 
to deduce 7, as others [69] have done. For the cold gases the 
precise value of 7, and its T-dependence are not particularly 
important, as long as it is non-zero at finite T. In this paper 
we will deduce reasonable values for 7(T) and Eo from to- 
mographic RF experiments. 



A. Extension Above Tc 

We can expand tpg{Q) at small Q, and in the normal state 
to find 



C(0) = ZMp„r = C/"'+x(0) 



(16) 



where the residue Z and pair dispersion (not shown) fig, are 
then determined ll76[l . This is associated with the normal state 
gap equation 

U-'+J2'-^^jr^=Zl^pa.r, (17) 
k 

Similarly, above Tc, the pseudogap contribution to A^ (T) = 

A2^(r) + A2 (T) is given by 



(18) 
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The number equation remains unchanged. In summary, when 
the temperature is above T^, the order parameter is zero, and 
A = Apg. Since there is no condensate, fipair is nonzero, and 
the gap equation is modified. From these equations, one can 
determine /x, A and jipair- 

B. Incorrect Criticism from the Drummond Group 

The Drummond group has made a number of incor- 
rect claims about our past work which we address here. The 
authors claim to have numerically studied the behavior associ- 
ated with the three possible pair susceptibilities. We note that 
there is no elemental numerical data in their paper, nor do they 
present details beyond their use of an "adaptive step Fourier 
transform" algorithm. This should be compared with work 
by the Tremblay group eHHIItI] and others |36, 57]. It is 
hoped that in future they will present plots of the t-matrix and 
self energy to the community, to the same degree that we have 
shared the output of our numerical schemes in References ifssll 
and 1 56]. Important will be their counterparts to Figs. 8a and 
9 (lower inset) in Ref. [55J, which show how rehable the form 
in Eq. ( fT2l i is for the full GGq self energy. More specifically: 
they have argued that the "decomposition into pg and sc con- 
tributions [see Eq. (fTTl l above], omits important features of 
the full theory". This claim is incorrect and is based on their 
Fig. 1 of Ref. 1 15] which can be seen to be unrelated to the 
pg and sc decomposition, since their analysis of our so-called 
"pseudogap theory" is confined to the normal phase. The de- 
composition only applies below Tc- Moreover, we refute the 
argument that the decomposition into sc and pg terms shown 
in Eq. ( fTTT l above is unphysical. This decomposition is asso- 
ciated with the fact that there are necessarily both condensed 
and non-condensed pairs in the Fermi gases at unitarity. This 
break-up is standard in studies of Bose gases. The details of 
how to describe the pg contribution, but not its necessary pres- 
ence in a decomposed fashion, are what varies from theory to 
theory. 

Importantly, the "discrepancies" associated with thermody- 
namical plots based on our approach should be attributed to 
the absence of a Hartree term, not to any deeper physics. The 
reader can see that if the usual (3 parameter is changed from 
—0.41 to around —0.6 the BCS-Leggett curve will be aligned 
with the others. 



IV. COMPARISONS OF THE SPECTRAL FUNCTION 

The fermionic spectral function is given by A(ti>,k) = 
—2lm[G , k)], where Gu is the retarded Green's function. 
In this section we want to explore its behavior both as a func- 
tion of wavevector k and of frequency uj. Of particular in- 
terest is the question of whether there is a pseudogap in the 
spectral function. There are different criteria for arriving at an 
answer. Importantly, depending of this choice the answer will 
be different for NSR theory (and also, it appears for the FLEX 
theory of Ref. 1 16]). The following definitions come from dif- 
ferent measurements of HTSCs which are not internally con- 




Figure 2: (Color online) Spectral function obtained from NSR theory 
(left column) and from the extended BCS-Leggett theory (right col- 
umn) at unitarity. Temperatures (T/Tf) from top to bottom are: (a) 
and (b) 0.24, (c) and (d) 0.34, (e) and (f) 0.55. The ranges of k/kp 
and uj/Ef are (0, 2) and (—2, 3), respectively. 



tradictory. Following Ref. |78] and references therein, we ex- 
amine two criteria for the pseudogap in HTSCs. 

1 . One can define the existence of a pseudogap as associ- 
ated with the observation that A(k, uj) as a function of 
UJ at k = kp exhibits a two-peak structure in the normal 
state. This definition is particularly useful for spectro- 
scopies such as ARPES which can probe the spectral 
function near the Fermi surface. 

2. Alternatively, the existence of a pseudogap can be iden- 
tified when the density of states (DOS) (which repre- 
sents an integral over k of the spectral function) is de- 
pleted near the Fermi energy. This definition appeals to 
tunneling experiments where the DOS can be measured. 

In addressing these criteria it is useful to refer to the 
spectral function of the BCS-Leggett ground state given by 




(1/2) [1 ± [S^k/Ek)]- As a function of frequency, there are two 
branches: the upper branch located at uj — Ei^ has weight 
ul and the lower branch located at oj = —Ek has weight w^. 
Since > A, the spectral function is gapped at all k. One 
recognizes two features in Abcs- First, there is particle-hole 
mixing which results in the two branches. Second, there is an 
upwardly dispersing and a downwardly dispersing symmetric 
contribution to the spectral function arising from the ± signs 
in Eq. Q. This is symmetric about the non-interacting Fermi 
energy. At finite temperatures, as one will see, both NSR the- 
ory and the extended-BCS-Leggett theory show particle-hole 
mixing in the sense that there are two branches in the spectra. 
In contrast, the fermionic dispersion in NSR theory does not 
lead to two symmetric upwardly and downwardly dispersing 
branches. The behavior of the finite T spectral function asso- 
ciated with BCS-Leggett theory, given in Eq. (|2]i is, however, 
rather similar to its superfluid analogue. 

It is unlikely that Eq. (|2]l will be appropriate at sufficiently 
high temperatures. Indeed, one can see from Figure 3 in Ref- 
erence [56] and the surrounding discussion that numerical cal- 
culations show this approximation is appropriate up to some 
temperatures of the order of T/Tc « 1.3 for a system near 
unitarity and in the absence of a trap. In a fully consistent nu- 
merical calculation one expects that as T is raised the pseudo- 
gap will decrease so that the pair susceptibility of the extended 
BCS-Leggett theory should eventually evolve from GGq to 
GqGq. In this way the fully numerical NSR scheme should 
be very reasonable at sufficiently high T, where the pseudo- 
gap begins to break down. Physically we can argue that the 
BCS-Leggett scheme is better suited to treating pairs which 
have pre-dominantly low momentum, and thus, it should ap- 
ply closer to condensation. For the purposes of comparison, 
in this section we apply Eq. (|2]i up to somewhat higher tem- 
peratures than appears strictly feasible. 

A very important physical distinction emerges between the 
different models for the pair susceptibility which is then re- 
flected in the fermionic self energy and ultimately in the spec- 
tral function. Because a dressed Green's function appears in 
BCS-Leggett theory, the t-matrix t{Q) at small q has a notably 
different behavior, particularly at low ui as compared with the 
NSR case. This is seen most clearly by comparing Figure 
2 and Figure 9 in Ref. llssll . This difference can be seen as 
a gap in the GGq t-matrix which serves to stabilize the pair 
excitations. In the normal state the pairs live longer when a 
pseudogap is present because of this feedback. As a result 
the behavior of the fermionic self energy is different, leading 
to a reasonable fit to Eq. (|2]i in GGq theory as shown by the 
lower inset in Figure 9 of Reference llssll . as compared with 
the poorer fit to Eq. (|2|i found in NSR theory and shown in 
Figure 8a from Ref. [55J. We will reach qualitatively similar 
conclusions in the next section of the paper. We summarize by 
noting that the extended BCS-Leggett theory focuses on low 
q pairs which dominate near condensation. NSR theory treats 
pairing without singling out low q only. Each of these the- 
ories should be appropriate in different temperature regimes 
of the normal state. Concomitantly, because of the enhanced 
stability of the pairs, the broadening of the spectral peaks will 
be considerably smaller in BCS-Leggett theory as compared 



k/kp 

Figure 3: (Color online) The plot of the function 
k'^f{ui)A{k.,oj)/2TT^ calculated from NSR theory (left col- 
umn) and from the extended BCS-Leggett theory (right column) at 
l/kpa = 0. Temperatures (T/Tf) from top to bottom are: (a) and 
(b) 0.24, (c) and (d) 0.34, (e) and (f) 0.55. The ranges of k/kp and 
uo/Ef are (0, 2) and (—3, 3), respectively. 



with NSR theory. 

To make a connection with experiments on ultra-cold Fermi 
gases, we regularize the attractive coupling constant via 
U^^ = m/{Ana) — X]k("^/^^)- choose as our units the 
Fermi energy Ep, or, as appropriate the Fermi temperature 
Tp, or Fermi momentum kp of a non-interacting Fermi gas 
with the same particle density. The unitary point where a di- 
verges is of particular interest because two-body bound states 
emerge. Since many-body effects renormalize the coupling 
constant, the fermion chemical potential remains positive at 
unitarity in both NSR theory and the BCS-Leggett theory. 
This implies that bound states in a many-body sense have not 
fully emerged. In our numerics, we choose 7(T) to be very 
roughly consistent with RF experiments. For the unpolarized 
case we set ^jEp — 0.12(T/rc) at unitarity and included a 
small background imaginary term Yj^/Ep — 0.05. 
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A. Comparison of Spectral functions via Contour plots 

Figure |2] present a plot of the spectral function at unitar- 
ity (l/kpo, = 0) obtained from NSR theory (left column) 
and from the BCS-Leggett t-matrix theory (right column) 
at selected temperatures. In the BCS-Leggett case we use 
the approximation that the self energy associated with non- 
condensed pairs is of a broadened BCS form (in Eq. (|2]i). The 
transition temperatures Tc/Tp = 0.238 and Tc/Tp = 0.26 
are obtained for the NSR and BCS-Leggett cases respectively. 
Both theories yield higher values than found 1 10] in quan- 
tum Monte Carlo simulations, where Tc/Tp ~ 0.15 at uni- 
tarity. The Tc curves in BCS-BEC crossover of NSR and the 
extended BCS-Leggett theories are shown in Ref. 1 14]. 

In Fig. 12] the comparisons are made at three different tem- 
peratures. The horizontal and vertical axes on each panel cor- 
respond to wave number and frequency and what is plotted 
in the contour plots is the fermionic spectral function for a 
three-dimensional homogeneous gas. The white areas cor- 
respond to peaks in the spectral function and they map out 
a dispersion for the fermionic excitations. With the possi- 
ble exception of the highest T NSR case (lower left figure) 
the spectral functions in all cases shown in Fig.|2]are gapped 
at small k/kp, which indicates the existence of particle-hole 
mixing The lower branch of the spectral function from the 
BCS-Leggett t-matrix theory clearly shows a downward bend- 
ing for k > kp which is associated with a broadened BCS-like 
behavior The spectral function of a phenomenological pseu- 
dogap model presented in Ref. ]65], which can be derived mi- 
croscopically from the BCS-Leggett approach, exhibits similar 
contour plots as ours from the extended BCS-Leggett theory. 
By contrast, in NSR theory the lower branch corresponds to 
very broad and very small peaks when k > kp which are 
barely observable. We see no clear evidence of a downward 
dispersing branch even at the lowest temperature above Tc- 

As a function of temperature, in the NSR case, the physics 
suggests a smooth evolution with increasing T towards a sin- 
gle branch, upwards dispersing- almost Fermi liquid disper- 
sion. It seems likely that as T is raised and the pairing gap 
becomes less important the pair susceptibility of the BCS- 
Leggett state should cross from GGq to GqGq so that the two 
schemes merge. This means that our previous simplification 
of the fermionic self energy S (as a broadened BCS form, see 
Eq. ( fT2] i) is no longer suitable in this high T regime. There are 
not really two different types of pseudogap, but rather the ex- 
tended BCS-Leggett theory represents pairs which are close to 
condensation- and thus predominantly low momentum pairs. 
[This is built into the approximation Eq. ( fT2l i. which was used 
to model the pseudogap self energy.] By contrast the NSR 
case considers pairs with a broad range of momenta. 

One can see that the vanishing of the pseudogap as temper- 
ature increases is different in the two theories. In the extended 
BCS-Leggett theory, the two branches approach each other 
and the gap closes while in NSR theory, the spectral func- 
tion fills in the gapped regime and its overall shape evolves 
toward a single parabola. We note that our spectral func- 
tion from NSR theory at unitarity looks identical to the re- 
sults in Ref. 11 2n] . even though the latter was computed with a 



more self consistent number equation. This implies that the 
difference between the two number equations (Eq. (|5]) and 
n = 2 G{K)) has no qualitative impact. Interestingly, 
the spectral function at the highest T from NSR theory resem- 
bles that presented in Ref. 1,161 for the near-Tc normal phase 
of GG t-matrix theory. 

It should be noted that this spectral function is not the 
quantity directly measured in momentum resolved RF stud- 
ies ll64ll . Rather what is measured there is the function 
/c^/(a;)A(k, w)/27r'^ which is plotted in Figure[3] This convo- 
lution preserves the large momentum part of the lower branch 
of the spectral function and suppresses the remainder. The 
downward bending behavior is clearly observed in the ex- 
tended BCS-Leggett theory again (right column) for all three 
temperatures. In contrast, only the lowest temperature plot of 
NSR theory {T/Tp = 0.24, sHghtly above Tc) shows a weak 
downward bending at large momentum. This downward dis- 
persion, however, cannot be observed at higher temperatures 
(left column). In the actual momentum-resolved RF experi- 
ments of trapped Fermi gases ll64ll trap averages enter so that 
the actual curves are substantially broader fTV\. A clear sig- 
nature of downward dispersion in experimental data will help 
determine whether the pseudogap phase with noncondensed 
pairs behaves in a way which is similar to the HTSCs, where 
this feature has been reported ||68|] . It should be noted that it 
is this signature which has been used in Ref. f64ll to arrive at 
an indication of the presence of pairing. 

The issue of what constitutes the proper definition of a 
pseudogap is an important one and we turn next to the first 
definition we introduced above in which one requires that 
the spectral function A{Lo,'k) as a function of lu exhibits two 
peaks around k k kp. It is clearly seen that in the BCS 
Leggett approach this definition is met for all three curves ex- 
hibited in Figure |2l Because it is more difficult to establish 
this for the NSR results plotted in Fig.|2] in Fig.|4]we address 
this question more directly for two different temperatures and 
a range of k values near kp.li should be clear that at the lower 
temperature (which is slightly above Tc), a pseudogap is seen 
in NSR theory, although we have seen that the peak dispersion 
is not well described by Eq. Q. This pseudogap should not be 
viewed as a broadened BCS like feature. At the higher tem- 
perature shown by the bottom panel of Fig. |4] there appears to 
be no indication of a pseudogap according to the first defini- 
tion. Only a single peak is found in the spectral function near 
kp. 

We next explore analogous curves in the BEC regime and 
thereby investigate how many-body bound states affect the 
distribution of weight in the spectral function. Figure |5] il- 
lustrates the spectral function on the BEC side of resonance 
with 1/kpa = 1 obtained from NSR theory (left column) and 
from the extended BCS-Leggett theory (right column) at se- 
lected temperatures. Here Tc/Tp — 0.22 for the NSR case 
andTc/Tp = 0.21 in the BCS-Leggett scheme. There is noise 
at small k in the spectral function of NSR theory which is pre- 
sumably a numerical artifact. In the presence of bound states 
(in a many -body sense), the lower branch of the spectral func- 
tion of NSR theory shows a downward bending near Tc, but 
it can be seen that this behavior rapidly evolves to an upward 
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co/Ep 

Figure 4: Behavior of the frequency dependent spectral function at 
unitarity in NSR case for (a) T/Tf = 0.24 and (b) T /Tp = 0.34 
for various wave-vectors k. This Figure suggests that the two-peak 
structure near kp associated with this crossover theory barely meets 
the (most restrictive) definition for the presence of a pseudogap near 
Tc. Away from Tc, the two-peak structure near is virtually not 
observable. 



dispersion as T increases (Fig.|5je)). In contrast, the spectral 
function from the extended BCS-Leggett theory exhibits the 
downward dispersion at all T indicated, rather similar to the 
behavior in the superfluid phase. 

One can further see from the figure that in the extended 
BCS-Leggett theory the lower branch has a much weaker 
spectral weight compared to that of the upper branch. This 
derives from the same phenomenon as in the BCS-Leggett 
ground state, where the coefficient becomes negligibly 
small in the BEC regime. The behavior of the NSR spec- 
tral function is rather different from its counterpart in BCS- 
Leggett theory at all three temperatures. If one plots the spec- 
tral function at fc = fcj? as a function of oj two peaks are 
present with the upper peak much sharper and narrower than 
the lower, as also reported in Ref . I,73ll . 




k/kp 



Figure 5: (Color online) Spectral function obtained from NSR the- 
ory (left column) and from the extended BCS-Leggett theory (right 
column) at 1/kpa = 1. Temperatures {T /Tp) from top to bottom 
are: (a) and (b) 0.25, (c) and (d) 0.34, (e) and (f) 0.55. The ranges of 
k/kp and uo/Ep are (0, 2) and (—4, 4), respectively. 

B. Comparison of Density of States 

We turn now to the DOS which, when depleted around the 
Fermi energy, provides a second criterion for the existence of 
the pseudogap. The DOS is given by 

iVH = ^^(c^,k). (19) 

k 

In the HTSCs, an above-Tc- depletion in the DOS around the 
Fermi energy, measured in tunneling experiments, provided a 
clear signature of the pseudogap ([78J and references therein). 
By contrast, in the ultra-cold Fermi gases the DOS has not 
been directly measured, although it is useful to discuss it here 
in an abstract sense. 

Figure |6] shows the DOS from the two theories at unitarity 
and at selected temperatures. The DOS based on the extended 
BCS-Leggett theory clearly shows a pseudogap at all three 
selected temperatures. Similarly, the DOS from NSR theory 
show a clear depletion around the Fermi energy {uj — 0) at 
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T/T^O.24 
T/T^O.34 



^ 




Figure 6: (Color online) DOS at unitarity from (a) NSR theory and 
(b) the extended BCS-Leggett theory. (Black) solid line, (red) dashed 
line, and (green) dot-dash line correspond to T /Tf ~ 0.24, 0.34, 
and 0.55. 



T/Tf = 0.24. This depletion is barely visible at 0.34. At 
higher temperature (T/Tp = 0.55), the depletion does not 
appear, and one sees only an asymmetric background. We 
note that our NSR results are similar to those in Ref. iTSil al- 
though different number equations were employed. With this 
criterion for a pseudogap one would conclude that NSR theory 
does have a pseudogap — at least at T « Tc- It is somewhat 
unlikely that the FLEX scheme (which at seems to behave 
similarly to the highest T NSR figures) a pseudogap would be 
present — via this second definition. 



C. Comparison of RF Spectra of Unpolarized Fermi Gases 

The RF current at detuning v also depends on an integral in- 
volving the fermionic spectral function. The current obtained 
from linear response theory is given by 1,541 



27r 
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where, for the present purposes we ignore the complications 
from final-state effects |J9i,,801 . as would be reasonable for the 



Figure 7: (Color online) RF spectrum at l/kpa = 0. The (black) 
solid dots, (red) solid lines, and (green) dashed lines correspond to 
the RF currents obtained from the experimental data (|61], supple- 
mental materials), the extended BCS-Leggett theory, and NSR the- 
ory. The values of T/Tf are (a)0.2, (b)0.22 (0.24 for the curve from 
NSR theory), (c)0.34, and (d)0.55. 



so-called "13" superfluidof ^Li. Here jTfcl^ is a tunneling ma- 
trix element which is taken to be a constant. The data points 
in Figure |2]correspond to measured tomographic spectra from 
Ref. II6II1 in units of the local Fermi energy, (see the Supple- 
mental Materials). The results from NSR theory are indicated 
by the dashed lines and from the extended BCS-Leggett the- 



12 



ory (GGq) by the solid (red) lines. To compare and contrast 
the spectra, we normalize each curve by its maximum and 
align the maxima on the horizontal axis, so as to effectively 
include Hartree shifts. The experimental data were taken at 
T/Tp = 0.2, 0.22, 0.34, 0.55. The RF spectra from the ex- 
tended BCS-Leggett theory are calculated at the same set of 
temperatures. The RF spectra from NSR theory, which are 
restricted to the normal phase correspond to the three higher 
temperatures: T/Tp = 0.24, 0.34, 0.55. 

The RF spectra from the extended BCS-Leggett theory at 
high temperatures indicate a double-peak structure, which was 
addressed in Ref. 1180(1 . This peak at negative RF detuning 
emerges at finite temperatures in BCS-Leggett theory as a re- 
sult of thermally excited quasiparticles. With increasing T, 
the weight under this peak increases although the peak-to- 
peak separation will decrease, following the temperature de- 
pendent pairing gap, as seen in the figure. When temperature 
increases, the peak at negative RF detuning grows and nearly 
merges with the peak at positive RF detuning so that it may 
not be resolved experimentally. 

By contrast, the RF spectra in NSR theory show a single 
peak which is broader than the experimental RF spectra. This 
is to be expected based on our analysis of the NSR spectral 
function in the previous section, where we saw that the sym- 
metrical upward and downward dispersing branches of BCS 
theory were not present. The RF spectra presented in Ref. B16ll 
using GG t-matrix (FLEX) theory also shows a broader (than 
experiment) single peak. 

In view of the contrast between the BCS-Leggett curves 
and experiment, it is natural to ask why is there no indica- 
tion of the negative detuning peak in these unpolarized exper- 
iments? One can contemplate whether this stems from the fact 
that (owing to large 7) the two peaks simply aren't resolved. 
This would yield a figure closer to that obtained from NSR- 
based calculations, which is associated with a rather broad 
peak structure. As a result it would not lead to a more satisfy- 
ing fit to experiment. At this stage we have no clear answer, 
but it will be important to investigate, as we do below, very 
slightly polarized gases to gain some insight. 

In Figure m we present the comparison on the BEC side of 
resonance. Here 1/kpa — 1 and T/Tp — 0.34. In this case 
both spectra show a double-peak structure. For convenience, 
we have scaled both spectra to their maxima and aligned the 
maxima. For the NSR case, the double-peaked feature reflects 
the negative fermionic chemical potential which is associated 
with bound states. It similarly reflects the stronger spectral 
weight of the upper branch in the spectral function which can 
also be examined in Figure |5] Our RF spectra in NSR theory 
are consistent with those presented in Ref. t70l] . 



V. EXTENDED BCS-LEGGETT THEORY: SIGNATURE OF 
Tc IN RF SPECTRUM 

We have tried in the paper to emphasize comparisons when- 
ever possible, but there are instances where other crossover 
theories (besides that based on the BCS-Leggett theory) have 
no counterpart. In the first of these we investigate the signa- 




v/£. 



Figure 8: (Color online) RF spectrum at l/kpa = 1. The (black) 
dot-dash line and (red) dashed line correspond to the RF currents 
obtained from the extended BCS-Leggett {GGq) theory and NSR 
theory at T/Tf =0.34. 



ture of the second order transition which should be a subtle, 
but nevertheless thermodynamically required feature of any 
crossover theory. The experimental RF spectra in Ref. ll6Tll62ll 
imply that the RF spectrum is more sensitive to the exis- 
tence of pairing rather than to superfluidity. That it evolves 
smoothly across Tc is due to the presence of noncondensed 
pairs. The extended BCS-Leggett theory has the important 
advantage in that it describes a smooth transition across Tc 
and should be a suitable theory for investigating this question. 
In contrast, NSR theory and its generalization below Tc, as 
well as the FLEX or Luttinger-Ward approach lfl6ll encounter 
unphysical discontinuities 

In the following we search for signatures of Tc in the RF 
spectrum as obtained from BCS-Leggett theory near Tc- Here, 
in contrast to Ref. (6f], we use constraints provided by our 
semi-quantitative fits to RF spectra (associated with the esti- 
mated size of 7(T)) to obtain a more direct assessment of how 
important these superfluid signatures should be. 

Figure |9] presents a plot suggesting how one might expect 
to see signatures of coherence in a tomographic (but momen- 
tum integrated) RF probe, such as pioneered by the MIT group 
ll6lll . It shows the RF current versus temperature at three dif- 
ferent detuning frequencies. The inset plots the RF character- 
istics indicating where the frequencies are chosen. One can 
see that there is a feature at Tc, as expected. This shows up 
more clearly in the lower figure which plots the temperature 
derivative. The same sort of feature has to be contained in the 
specific heat Isill which represents an integral over the spec- 
tral function. What does it come from, since RF is not a phase 
sensitive probe? The feature comes from the presence of a 
condensate below Tc- What is distinguishing condensed from 
non-condensed pairs is their self energy contribution. In the 
HTSCs |74] and also in the BCS-Leggett formulation the self 
energy from the non-condensed pairs is taken to be of a broad- 
ened BCS form in Eq. ( fT2] i. By contrast, the non-condensed 
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Figure 9: (Color online) (a) RF current as a function of T at 
l/kpa — for v\/Ef = 0.15 (black dashed line) and v\/Ef ~ 
0.1 (red solid line) obtained from the extended BCS-Leggett theory. 
Inset: RF currents as a function of detuning at T/Tf ~ 0.22 (dot- 
dash line), 0.26 (dashed line), 0.3 (solid line). The two arrows in- 
dicate and U2. (b) The slopes of the RF currents from (a). The 
vertical dashed line indicates Tc- 



pairs live infinitely long and so have no damping 7. These 
are the effects which are represented in Figure |9] In this way 
the figure shows that there are features at Tc which can in 
principle help to distinguish the ordered state from the normal 
pseudogap phase. 



VI. EXTENDED BCS-LEGGETT THEORY: RF 
SPECTRUM OF POLARIZED FERMI GASES 



Another strength of the BCS-Leggett approach is that it can 
address polarized gases at unitarity, which are not as readily 
treated ll40l l4lll in the alternative crossover theories. In Fig- 
ure[TO]we plot the RF spectra from the extended BCS-Leggett 
theory and the experimental RF spectra from Ref. [61]. Since 
the experimental RF spectra were obtained from RF tomog- 
raphy of trapped polarized Fermi gases, we follow a similar 
procedure to extract our RF spectra at varying, but comparable 
locations from a similar trap profile. Also indicated are the po- 
larizations p. If we make fewer restrictions on the choice of ra- 
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Figure 10: (Color online) RF spectra of polarized Fermi gases 
at unitarity. (Black) dots and (red) solid lines correspond to the 
experimental data from Ref. lIcT] and results from the extended 
BCS-Leggett theory. The left (right) column shows the RF spec- 
tra for the majority (minority) species. The local temperature 
T/Tfi and local polarization p for the experimental data (ex) 
and the theoretical results (th), {Tex /Tf-\ ,Tth/TF^ ,Pex,Pth), are: 
(0.05, 0.04, -0.04, 0) for (a) and (d); (0.06, 0.13, 0.03, 0.07) for (b) 
and (e); (0.06, 0.15, 0.19, 0.23) for (c) and (f). 



dial variable, the agreement is better as is shown in Ref. 162 
To compare the results, we normalized the maxima and align 
the spectra, thereby introducing a fit to Hartree contributions. 
The left (right) column shows the RF spectra of the majority 
(minority). Here we set -i/Ep = 0.05 and Eq/Ef = 0.02. 

The experimental data points from the left hand column can 
be compared with those in Figure [7] which are for the p = 
case, and it is seen that even at very small polarizations (say 
p « 0.03) the negative detuning peak becomes visible. In- 
deed, it appears here to be larger than the theoretically esti- 
mated negative detuning peak height. A possible explanation 
for why the double-peak structure can be resolved experimen- 
tally in polarized but not in unpolarized Fermi gases is because 
the existence of excess majority fermions causes a negative 
RF-detuning peak even at low temperatures. At these lower T 
the separation between the two peaks can be large in the exper- 
imental RF spectra of polarized Fermi gases. In contrast, for 
an unpolarized Fermi gas the negative RF-detuning peak due 
to thermally excited quasiparticles only becomes significant at 
high temperatures around Tc and above. Here the separation 
between the two peaks may not be as readily resolved. As ex- 
pected, at low temperatures there is only a single peak in the 
RF spectra of the minority. We notice that at extremely high 
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Figure 1 1 : (Color online) Reproduction of Fig.5 in the Supplemental 
Materials of Ref. |72]. Red curves correspond to experimental data 
from Fig. \T0\ (a) and (b). Black curves are RF currents calculated 
from BCS-Leggett theory. 



polarization, polaron-like behavior has been observed in RF 
experiments 1,82,1 . whose explanation has attracted a great deal 
of attention in the theoretical community Issilsill . These ef- 
fects have not been incorporated in our BCS-Leggett formal- 
ism, where the normal state has been assumed to be strictly 
non-interacting jssll . 

The Ketterle group [l72t] has argued that it should be possi- 
ble to extract the pairing gap size from RF spectroscopy in po- 
larized gases at very low temperatures. In Fig.[TT]we present a 
plot from their paper (Supplementary Materials) which relates 
to their procedure. This figure presents a fit to a generalized 
BCS-Leggett ground state in the presence of polarization. The 
red curves correspond to the actual data and the black curves 
are obtained from this theory. An additional resolution broad- 
ening is included in the theory and one can see that this theo- 
retical approach appears to be in quite reasonable agreement 
with experiment. In this way there is some support for this 
simplest of ground states — at least in the polarized case. 



VII. CONCLUSION 

The goal of this paper is to communicate that BCS-BEC 
crossover theories are very exciting. They are currently being 
clarified and developed hand in hand with experiment. For 
the Fermi superfluids, unlike their Bose counterparts, we have 
no ready-made theory. In this paper we confine our attention 
to the normal phase, although we have presented a discus- 
sion of some of the controversial issues which have surface 
in the literature below Tc- We view the principal value of 
this paper is the presentation of comparisons of two differ- 
ent crossover theories and the identification of (mostly future) 
experiments which can help distinguish them. The two the- 
ories we consider are the extended BCS-Leggett theory and 
that of Nozieres and Schmitt-Rink. We chose not to dis- 
cuss the FLEX or Luttinger-Ward scheme in any detail be- 
cause it is discussed elsewhere 1.161 , and because there are 
concerns that, by ignoring vertex corrections, this approach 
has omitted the important physics associated with the pseu- 
dogap. These concerns are longstanding ll22l |24, 35, 36]. 
Here we have argued that the extended BCS-Leggett the- 
ory is the one theory which preserves (broadened) BCS fea- 
tures into the normal state over a significant range of tem- 
peratures. Even above Tc one finds that the fermionic exci- 
tations have an (albeit, smeared out) dispersion of the form 

i?k « ±^{ek " m)^ + in the normal state. We find that 
NSR theory does not have this dispersion, although it has a 
pseudogap by all other measures. Interestingly high Tc super- 
conductors have been shown to have this dispersion in thek 
normal state [68] and it is generally believed i69[|74ll75ll that 
their fermionic self energy can be fit to a broadened (c?-wave) 
BCS form T.pg{K) « Apg(k)V(w + ek - M + h)- 

In this paper we show that one can identify both physically 
and mathematically the difference between the two normal 
states of the different crossover theories. Mathematically be- 
cause BCS theory involves one dressed Green's function in 
the pair susceptibility, it leads to a low frequency gap in the 
t-matrix or pair propagator (at low q). Physically this serves 
to stabilize low momentum pairs. This helps us to understand 
that the pseudogap of NSR theory does not incorporate pri- 
marily low momentum pairs, but rather pairs of all momen- 
tum and that it should be better further from condensation. 
Indeed, this is reenforced by our observation that at higher T, 
feedback effects which distinguish the two theories becomes 
less and less important and the BCS-Leggett pair suscepti- 
bility, GGq, crosses over to something closer to GqGq as in 
NSR theory. Our simplest approximation for the self energy 
in Eq. (|2|i is no longer suitable once temperature exceeds, say 
T/Tc ~ 1.5. Indeed this is reenforced by earlier numerical 
observations ll55il56ll . 

As a result, we believe that both theories are right but in dif- 
ferent temperature regimes. Moreover, this serves to elucidate 
another concern about NSR theory (and FLEX theory)- that 
they are associated with an unphysical first order transition. 
Both theories change discontinuously in going from above to 
below Tc- In the superfluid phase the coupling which is in- 
cluded in all other theories is between the non-condensed pairs 
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and the collective modes of the condensate, even though in the 
normal state one couples the fermions and the non-condensed 
pairs. In the extended BCS-Leggett theory, (as seems reason- 
able, in the vicinity of Tc both above and below), the dom- 
inant coupling is, indeed, between non-condensed pairs and 
fermions. These (effectively, pseudogap) effects will behave 
smoothly across Tc- The Goldstone modes which turn on at 
Tc are highly damped in its vicinity, where the condensate is 
weak. Only at lower T should their coupling become the more 
important. 

In summary, a central conclusion of this study of the spec- 
tral functions of the extended BCS-Leggett theory and NSR 
theory is that one may expect that the former is suitable near 
Tc due to its similarity to BCS theory while the latter better 
describes the normal phase at much higher T as the system 
approaches a Fermi liquid, and concomitantly, the pseudogap 
begins to disappear. In the course of this work we have found 



that the theoretical RF spectra from both theories agree (only 
semi-quantitatively) to about the same extent with experimen- 
tal data at unitarity. The BCS-Leggett approach has the advan- 
tage that it can address the RF spectrum of generally polarized 
Fermi gases without the problems which have been noted [4 ill 
for the NSR approach. However, momentum resolved experi- 
ments may be the ultimate way of distinguishing experi- 
mentally between different theories. 
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